Quantum dynamics of a single, mobile spin impurity 
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Quantum magnetism describes the properties of many materials such as transition metal oxides 
and cuprate superconductors. One of its elementary processes is the propagation of spin excita- 
tions. Here we study the quantum dynamics of a deterministically created spin-impurity atom, as 
it propagates in a one-dimensional lattice system. We probe the full spatial probability distribution 
of the impurity at different times using single-site-resolved imaging of bosonic atoms in an optical 
lattice. In the Mott-insulating regime, a post-selection of the data allows to reduce the effect of 
temperature, giving access to a space- and time-resolved measurement of the quantum-coherent 
propagation of a magnetic excitation in the Heisenberg model. Extending the study to the bath's 
superfluid regime, we determine quantitatively how the bath strongly affects the motion of the im- 
purity. The experimental data shows a remarkable agreement with theoretical predictions allowing 
us to determine the effect of temperature on the coherence and velocity of impurity motion. Our 
results pave the way for a new approach to study quantum magnetism, mobile impurities in quantum 
fluids, and polarons in lattice systems. 



Deepening our knowledge of quantum magnetism [1] 
is one of the most important goals in quantum simu- 
lation. In particular, one highly desired aim is a bet- 
ter understanding of high-T c cuprate superconductors, 
which are believed to be described by Heisenberg-type 
effective spin models in the limit of low doping |2J. 
More generally, the physics characterized by the Heis- 
enberg model governs the properties of many strongly 
correlated materials such as transition metal oxides [3], 
and allows the realization of various remarkable types 
of spin order ranging from spin solids to spin liquids 
[TJ |4]. In low-dimensional quantum magnets, particu- 
larly rich physics emerges due to the dynamics of spin 
excitations |5J. A basic mechanism of quantum mag- 
netism in strongly correlated electronic systems is su- 
perexchange, in which opposite spins on adjacent lat- 
tice sites coherently exchange their positions. Ultracold 
atoms in optical lattices offer an ideal testbed to explore 
these phenomena in a controlled experimental environ- 
ment [6j, as shown by the observation of superexchange 
in double- well systems |7J or plaquettes |8J. Further- 
more, the recently demonstrated single- atom-resolved 
detection [9j [10] opens entirely new prospects for the 
quantum simulation of strongly correlated spin systems. 
For example, this technique has enabled the simulation of 
one-dimensional anti-ferromagnetic Ising spin chains via 
mapping of the site occupation onto a pseudo-spin 
In ultracold-atom experiments, however, observation of 
superexchange-based quantum magnetism in many-body 
systems has so far been hindered by the fact that typical 
temperatures are much larger than the superexchange 
coupling energy. 



In this work, we report on the first space- and time- 
resolved observation of a coherently propagating spin- 
wave, by tracking the motion of a deterministically cre- 
ated single-spin impurity in a one-dimensional (ID) spin 
chain (Fig. [I]). Specifically, we studied the quantum dy- 
namics of a mobile boson of type ' | ' on a ID lattice, sur- 
rounded by a bath of bosons of type ' *f '. Such a system 
can be described within a two-species single-band Bose- 
Hubbard model, parametrized by the spin-independent 
single-particle tunneling rate J and on-site interaction 
energy U (Supplementary Information). Deep in the 
Mott insulator (MI) regime (U ^> J) with unity filling, 
this model can be mapped to the isotropic Heisenberg 
model [IMS]: 
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in which the effective superexchange coupling J ex = 
4J 2 /U arises from a second-order tunneling process. The 
bosons ' t ' and ' I ' are identified with spin states | t) and 
| 1) and the corresponding operators of the effective spin 
system are S = (%, Sj, SJ), S+ = = | t);U \j , 
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with the number operators h a j for bosons of type a = t 
, I on lattice site j. In the case of a single flipped spin in a 
ferromagnetic domain, the second term of Eq. ( lb ), which 
describes the longitudinal spin coupling, only gives rise to 
an energy offset and can therefore be neglected. The re- 
maining first term of Eq. ( lb ) is structurally equivalent to 
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the single-species single-particle tunneling Hamiltonian 
within the tight binding model of a ID lattice system 
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Figure 1. Coherent propagation of a single spin excita- 
tion in the Heisenberg model, a, A single spin is nipped 
at the centre site of a ID spin chain. Each spin is coherently 
coupled to its neighbors via the superexchange coupling J ex - 
b, c Coherent evolution showing the polar angles of the spins 
after times 2ft/ J ex and Ah/ J ex - Circles visualize the quantum 
uncertainty of the transverse spin component. 



[TJ |5j [TB], with the tunneling rate J being replaced by 
J ex /2, and the atomic creation and annihilation operat- 
ors being replaced by the spin- flip operators Sj~ and S~ . 

For larger J/U, the mapping to the Heisenberg model 
breaks down due to the fluctuations of the site occupan- 
cies and due to higher-order processes, and no analytic 
solution describing the time evolution of the impurity ex- 
ists. In this regime, the motion of the impurity is mod- 
ified compared to both a free particle and an isolated 
magnetic excitation. This effect can be interpreted as a 
result of interactions between the impurity and low-lying 
density excitations, which resemble phonons, leading to 
polaron-like effects. 

In order to create ID spin chains in the optical lattice, 
we followed the experimental procedure of our previous 
work [17J. We started by preparing a two-dimensional de- 
generate gas of about 170 87 Rb atoms in a single antinode 
of a vertical optical lattice (lattice spacing a\ at = 532 nm) 
along the z direction. We then switched on two hori- 
zontal optical lattice beams within 120 ms with potential 
depths V x = 10.0(3) E r and V y = 30(2) E r (the num- 
ber in parenthesis denotes the uncertainty of the last 
digit), where E r = h 2 /(8maf at ) denotes the recoil en- 
ergy and m is the atomic mass of 8T Rb. This creates 
an array of parallel ID Mis along the x direction, con- 
taining each 8 to 16 atoms. The atoms were initially 
prepared in the hyperfine state \F = l,rap = —1) (= 
't')- We introduced the spin impurity by changing the 
hyperfine state of an atom at the centre of the chain to 
state \F = 2, = —2) (= ' I ') using single-site address- 
ing [18J. In this scheme, a a - -polarized, off-resonant laser 
beam at 787.65 nm wavelength focussed onto the selected 
lattice site results in a negative energy shift (attractive 
potential) only for the state ' 4- ' while leaving the initial 
spin state ' t ' almost unaffected. A microwave pulse, res- 
onant with the shifted atomic transition, then produced 
the spin- flip from 'f'to c |'. In contrast to our previous 



work [18J, a spatial light modulator generated an address- 
ing beam profile in form of a line instead of a circular- 
shaped Gaussian beam, in order to create the impurity 
in all ID chains simultaneously (Methods). This novel 
multiple-site addressing technique offers the advantage 
to prepare an arbitrary spin pattern in the lattice more 
rapidly. In addition, the line-shaped beam allowed us 
to hold the impurities in all ID systems, and to release 
them simultaneously by switching off the beam within 
1 ms. We then allowed the impurity in each chain to 
propagate during a variable hold time, and finally froze 
the dynamics by rapidly increasing the lattice depth of 
all axes to > 80 E r within 300 ps. 

The resulting spin distribution was detected by single- 
site-resolved fluorescence imaging using a high-resolution 
microscope objective [10]. We used two alternative detec- 
tion methods to determine the position of the impurity: 
(i) the direct imaging of the impurity spin component 
' 4 ' after removing all atoms in state ' t ' with a resonant 
laser pulse ('positive image') and (ii) detecting it as an 
empty site in the bath of ' t '-atoms after removing the 
spin impurity ('negative image'). The latter has the ad- 
vantage to provide information about the thermal excita- 
tions (holons and doublons) of the system as well, which 
are also detected as empty sites. Although we cannot 
distinguish between spin impurity and thermal excita- 
tions, post-selecting samples with only one empty site 
in the chain enables us to filter out a lower-temperature 
subset of the data. If the only empty site arises from 
a spin-flipped atom, the position of the impurity can be 
determined exactly and the selected samples contain no 
additional excitations. 

We first studied the time evolution of a spin impurity 
deep in the MI regime using negative images. The prob- 
ability distribution of its position was obtained by aver- 
aging data from different chains, post-selected to contain 
only one empty site within the central 10-14 sites. The 
distributions [see Fie;. H for J/U = 0.053(7)] show clear 
maxima and minima, resulting from the quantum inter- 
ference due to the coherent evolution of the spin impurity. 
Because of the weaker superexchange coupling, the ob- 
served dynamics occurs on timescales much longer than 
the tunneling time (h/ J = 4 ms) that would characterize 
the motion of non-interacting atoms. We compared our 
data to the time evolution of the spin impurity in the 
exactly solvable homogeneous Heisenberg model at zero 
temperature. In this model, the probability of finding the 
impurity at time t on site j after starting its evolution at 
t = from the centre of the chain (j = 0) is 
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where Jj is the Bessel function of the first kind [19J. 
A single fit to all distributions observed at different hold 
times (red curves in Fig.[2| with the superexchange coup- 
ling as a free parameter yields good agreement with the 
data for J ex /h = 65(1) Hz, which is close to J ex /^ = 
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Figure 2. Dynamics of a mobile spin impurity. The 

top panels in a and b show fluorescence images of the atoms 
(left) taken with the 'negative' imaging technique (see main 
text) together with the reconstructed atom distribution in the 
central region (right). The ID systems are oriented horizont- 
ally. The red vertical stripe denotes the initial position of 
the spin impurity (detected as an empty lattice site). For the 
generation of an effective low-temperature subset, only the 
samples containing exactly one empty site were kept (green 
tick marks) , while those containing more than one empty site 
were discarded (red crosses) . The histograms in a-d show the 
position distribution of the spin impurity after different hold 
times for J/U — 0.053. Each histogram is obtained from an 
average over 200-250 ID systems. The error bars denote the 
la statistical uncertainty. The red line is a simultaneous fit to 
all distributions with the analytic solution of the Heisenberg 
model of equation (J2|, yielding J ex /^ = 65(1) Hz. 



51 (1 8 ) Hz obtained from an ab-initio band-structure cal- 
culation using the independently measured lattice depths 
(see also Fig.[4|. Small differences between the data and 
the model result from the limited efficiency of the post- 
selection process (Supplementary Information). We note 
that the effect of the external potential on the spin dy- 
namics can be neglected in the MI limit, as long as the 
impurity has not yet reached the edge of the system (Sup- 
plementary Information) . 

To examine the effects of temperature, we measured 
the same distribution using positive imaging (Fig.[3|i), 
in which we directly observed the impurity position in 
the finite-temperature bath. Compared with the dis- 
tribution from the post-selected negative images with 
one empty site (Fig. [3^), it has almost the same width, 
but less contrast. This indicates that the propaga- 
tion velocity of the spin impurity is not much af- 
fected by temperature, although its coherent evolution 
is hindered. Our experimental observation is quantit- 
atively confirmed by numerical simulations (Fig. [3J3 and 
d) using the time-dependent density-matrix renormaliza- 
tion group (t-DMRG) algorithm, including temperature 



(Methods). Our intuitive interpretation of this effect is 
that thermal excitations can move back and forth across 
the spin impurity, introducing phase slips that wash out 
the coherence of the distribution. However, a single excit- 
ation alters the position of the spin impurity at most by 
one site in a ID system and as a consequence the width 
of the distribution is less affected. For a quantitative 
analysis, we evaluated the contrast and the propagation 
speed from the simulated distributions for different tem- 
peratures from T = to T = 0.2U/k B (Fig.|5]l). We 

defined the Contrast as C = (Pmax-^min)/ (^max+^min), 

where P m ax is the peak value at positions other than the 
centre, and P m i n denotes the minimum in the centre re- 
gion between the peaks (i.e., C = if there is no other 
maximum except for the central peak). The contrast ob- 
tained from the positive image equals that of the sim- 
ulated distributions at T = 0.15 U/k B (blue point in 
Fig. [3]i). This is consistent with the system's temper- 
ature of T = 0.14(3) U/k-B, determined independently 
by comparing the distribution of thermal excitations to 
an analytical model in the atomic limit [T0J. The con- 
trast from negative images corresponds toT = 0.11 U/k B 
(green point in Fig.|3]i), demonstrating that the post- 
selection enables us to extract a lower-temperature sub- 
set of the data. The actual temperature of this subset is, 
in fact, even lower, because our simulation does not take 
into account the finite spin- flip fidelity of 88(5)% and the 
fact that the initial position of the spin impurity could 
be empty due to a thermal excitation (see Supplement- 
ary Information). We also determined the propagation 
speed of the spin impurity as a function of temperature 
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Figure 3. Effect of thermal excitations on the coher- 
ent spin dynamics. Shown are position distributions of the 
spin impurity for J/U = 0.053 and a hold time of 40 ms, 
a generated from positive images (blue bars) and b after 
post-selection from negative images containing one empty site 
in the chain (green bars), c, Position distributions from t- 
DMRG simulations for different temperatures, d, Contrast of 
the simulated distributions (open circles) and propagation ve- 
locity (gray filled circles) as a function of the temperature T. 
The blue and green filled circles show the contrast extracted 
from the experimental data of a and b, respectively. Error 
bars denote the la uncertainty. 
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Figure 4. Spin dynamics across the superfluid-to-Mott-insulator transition, a-c, Spin impurity dynamics in the SF 
regime (J/U = 0.32) close to the critical point [(J/U) c w 0.3] for different hold times. The upper panels display fluorescence 
images of the impurity spins after removing the other spin component ('positive image'). The ID chains containing more than 
one atom were excluded from the data analysis. The white vertical stripe highlights the initial position of the nipped spin. The 
lower panels show the position distribution averaged over about 300 chains (blue bars) together with a t-DMRG simulation at 
T = 0.11 U/kB (red line), d, Measured velocities of the spin impurity for different values of J/U extracted from negative (green 
circles) and positive images (blue circles); horizontal and vertical error bars indicate the la uncertainties of the lattice depth 
and the combination of fit error and uncertainties of J, respectively. The dark gray line shows scaling with 4J 2 /[7, whereas 
the brown line indicates the propagation velocity of a single free particle (J/U = oo). The gray shaded region shows results 
from a t-DMRG simulation at T — taking into account varying initial atom numbers. The area denotes the la fit error to 
the simulated distributions. 



by fitting the simulated distribution with the function 
of Eq. (12), rewritten as Pj (t) = [Jj (^/ai at )] 2 , where we 
introduced v as the speed of the impurity. We found 
that the change in velocity within the simulated temper- 
ature range is only 4% (Fig.[3]i). This result is consistent 
with our observation of similar propagation speeds from 
positive and negative images (see also Fig.[4]i) and the 
intuitive argument mentioned above. 

In the MI regime, the impurity propagates in an en- 
vironment that is defect-free except for thermal excita- 
tions, leading to the coherent quantum motion described 
in the previous section. The situation is radically differ- 
ent in the superfluid (SF) regime, in which the bath of ' t ' 
spins contains many low-energy excitations as a result of 
quantum fluctuations of the on-site density. In one di- 
mension, those excitations are known to be phonon-like, 
which lead to a Frohlich-type Hamiltonian now describ- 
ing the drastic modification of the impurity motion in the 
bath by the propagation of a polaron [20J. In order to 
investigate this regime, we measured the local impurity 
density at different times for increasing values of J/U. As 
the impurity is always prepared inside a Mott insulating 
state at J/U = 0.053, we decreased the lattice depth V x 
along the chain within 50 ms to reach the desired final 
value of J/U. While preparing the bath in this way, we 
kept the impurity spin pinned by the addressing beam 
before releasing it. For large final J/U, close to the SF- 
MI critical point, post-selecting samples with only one 
empty site is no longer effective for tracking the spin be- 
cause of the proliferation of particle-hole pairs induced by 
quantum fluctuations [17J. We therefore recorded the dy- 
namics at large J/U using positive images only, in which 
we directly detected the impurity spins (Fig.[4^i-c). 



Since there is no full analytical solution describing the 
dynamics of the spin impurity for intermediate values 
of J/U, we chose to analyze the experimental data by 
fitting the width of the position distribution with the 
function Pj (t) with the velocity as an adjustable para- 
meter. This function is indeed the analytic solution 
both for the Heisenberg (MI) regime (J/U <C 1) with 
v = Jex^iat/fi- = 4(J/Z7) • (Jai a t/fi) and for a totally 
free impurity (J/U = oo) with v = 2Ja\ at /H. We found 
that a fit with Pj (t) captures well the edges of the posi- 
tion distribution, allowing us to determine the maximum 
propagation velocity for all J/U values used in this study. 

We used the same fit to determine the propagation ve- 
locity from the numerical simulations. These took into 
account the statistical atom number fluctuations from 
the experiment as well as the trapping potential, which 
are both important in the SF regime. Over the whole 
range of J/U accessed in the experiments, the measure- 
ment shows a good agreement with the simulation, both 
for the position distribution (Fig.^-c), as well as for the 
impurity velocity (Fig.[2]l). For small J/U, the experi- 
mentally determined spreading velocity normalized to the 
tunneling rate J increases as v/J = 4J/C7, as expected 
from the superexchange coupling J ex = 4 J 2 /U '. 

The spreading velocity becomes smaller than the ve- 
locity v/J = 4 J/U close to the point where the SF-MI 
transition would occur in a homogeneous system in the 
thermodynamic limit. We interpret this change of velo- 
city as a result from the interaction of the impurity with 
the increasingly large number of low energy excitations 
of the bath, which now affect not only the contrast but 
the velocity as well. In the SF regime, experiment and 
numerical simulations show that the normalized velocity 
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is only v/J^ a\ at /H, i.e., it is only about half the velocity 
expected for free-particle tunneling. This indicates that 
the system's strong interactions in the superfluid regime 
still slow down the motion of the impurity significantly. 
Such a behaviour is consistent with analytical calcula- 
tions [21 J and can be interpreted as a mass increase of 
the impurity which is expected from polaronic physics 
[HI22H25]. 

In conclusion, we have observed the space- and time- 
resolved dynamics of a deterministically created spin 
impurity in a ID lattice across the SF-MI transition. 
By post-selecting measurements we filtered out a lower- 
temperature subset of our data, allowing for the observa- 
tion of a coherently propagating magnetic excitation in 
the MI regime, consistent with magnetic superexchange. 
Across the MI-SF transition our measurement shows that 
the increasingly large number of low energy excitations 
of the bath affects both coherence and velocity of the 
impurity motion. The velocity, which can be related to 
the renormalized tunneling amplitude in the lattice, is 
significantly reduced compared to a free particle. This 
is consistent with the expected effect of such a bath and 
the mass increase due to polaronic effects, which is a field 
of intense ongoing research [22-29J. The good agree- 
ment with numerical results demonstrates that our ex- 
periment is suitable for the quantitative study of mobile 



impurities in Bose-Hubbard models [3D]. In particular, 
a predicted and fundamentally new universality class of 
single-particle excitations in ID systems could be tested 
by measurement of spin-flip response and Green's func- 
tion [21]. Other natural extensions of this work would be 
the investigation of the correlated propagation of several 
flipped spins due to the ferromagnetic spin-spin coupling 
in the Heisenberg model [31] or the quantum simulation 
of the impurity dynamics in 2D, where numerical simula- 
tions and analytical studies are much more difficult. Our 
multiple-site addressing technique could be used to pre- 
pare non-equilibrium states of quantum many-body sys- 
tems such as domain walls [32] and study their evolution. 
Further experiments could probe the fractionalization of 
excitations that naturally occurs in one dimension, either 
when magnons decay into two spinons, or when spin and 
charge excitations separate [20] . 
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METHODS 
Simultaneous multiple-site addressing 

We prepared the initial state with a single flipped spin 
at the centre of each one-dimensional (ID) system by us- 
ing a new technique to address several lattice sites sim- 
ultaneously. In our previous work, we used a single cir- 
cular shaped Gaussian laser beam focussed onto a single 
lattice site together with a microwave field for address- 
ing an individual spin in the lattice [15] . Instead of a 
circular beam, we now used arbitrary light-intensity pat- 
terns created with a spatial light modulator (SLM) of 
digital-micromirror-device (DMD) type in order to ad- 
dress several atoms at the same time. We generated a 
line-shaped light pattern by reflecting a Gaussian laser 
beam off the DMD (Texas Instruments, DLP Discovery 
4100, 1024x768 pixel, 13.7 /im pixel size), and coupled 
it into the high-resolution imaging setup with a dichroic 
mirror. The magnification was such that one pixel of 
the DMD corresponds to ai at /8 ~ 70 nm in the object 
plane. Compared to the diffraction limited spot size of 
~ 600 nm this oversampling allowed us to implement an 
error-diffusion algorithm [33] to generate beam profiles 
with about 10% root-mean-square variations of the peak 
intensity. Similar to our previous work [18J, we applied a 
feedback that shifted the pattern on the DMD for com- 
pensating slow phase drifts of the optical lattice. 

The addressing laser had a wavelength 787.65 nm and 
was <j~ polarized in order to minimize the light shift for 
atoms in the initial state \F = l,77ip = —1), while cre- 
ating a differential light shift between the initial state 
and \F = 2, Trip = —2) of ~ 35 kHz. In order to flip the 
spin, we used an HSl-pulse [18J of 20 ms duration, 40 kHz 
sweep width and —30 kHz frequency offset of the sweep 
centre from the bare resonance. The addressing beam 



was switched to full intensity for the spin-flip, and then 
back to 10% of the intensity for pinning the spin impurity 
(see main text) with s-shaped ramps of 10 ms duration. 
In relatively deep optical lattices (V x ,V y ,V z = 20 E r ) we 
obtained a spin- flip fidelity of 97{±l)%, which is the same 
as in our previous experiment [18J. In this work, the lat- 
tice depths were kept at a lower value along the ID chains 
(V x = 10 E r ,V y = 30E r ,V z = 20 E r ) in order to reduce 
heating from the lattice ramp after the spin preparation, 
and we achieved a spin- flip fidelity of 88(5)%. 



Numerical simulations 

In order to compute the time-dependent ID many- 
body state of the mobile impurity and the bath, we 
used numerical methods such as time-dependent matrix 
product states (t-MPS, a.k.a. time-dependent DMRG) 
[34H38] and simulated the two-species Bose-Hubbard 
model of Eq. (1) in the Supplementary Information. For 
our parameters, considering MPS of dimension 200 is suf- 
ficient to describe the impurity dynamics for J/U < 0.25, 
while for J/U > 0.25 a dimension of 400 is more ap- 
propriate, as the subsystem entanglement grows signific- 
antly close to and inside the SF regime. We incorporated 
the experimental conditions and preparation sequences 
as follows: The finite temperature of the initial MI state 
(experimentally prepared with an effective J/U = 0.053, 
at a lattice depth V x = 10E r: V y = 30 E r ) is considered 
by randomly sampling product states of localized atoms 
from the grand-canonical Gibbs-ensemble for the given 
initial U ', with J set to zero [10] . 

Deep in the MI regime, this is a very good approxim- 
ation. Chemical potential \i and external confinement 
<^trap were determined from the experiment [lOj . We 
generally used 100 randomly drawn initial states, which 
are then subjected to time evolution in accordance with 
the parameter changes in the experiment. To obtain 
the thermal average of the impurity dynamics, the time- 
dependent local densities for each initial state 
are used to form an unweighted average. 

When we simulated the system in the Heisenberg limit, 
deep in the MI regime, J/U does not change during the 
state preparation and the time evolution. The atom on 
the central site is spin-flipped and then evolved in time 
together with the background. When simulating the sys- 
tem at larger J/U than the initial one, we follow the 
experimental sequence exactly. Correctly modeling the 
release of the spin impurity is essential, as the times- 
cales for switching the addressing laser are similar to the 
tunneling times. Specifically, we simulated the adiabatic 
state preparation as it is done in the experiment using 
an initial ramp-down of the lattice depth V x within 50 ms 
according to 



V x (t)=Vi 



1 + coth 



tanh 



2t-T N 



> (3) 
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with Vi = 10 E r , T = 50 ms, r = 15 ms and V) the desired 
final lattice depth. We also incorporated into the time- 
dependent Hamiltonian parameters the switching off of 
the addressing laser within 1 ms. 



SUPPLEMENTARY INFORMATION 

A. The two-species Bose-Hubbard model and the 
Heisenberg model limit 

We described the effective dynamics in a ID lattice by a 
two-species single-band Bose-Hubbard model [5], which 
is valid at the temperatures, interaction strengths and 
lattice depths relevant to this study: 

#BH = - J ^ Qjb^k + ^2 ^nrj&r'J ~ ) 



(4) 



Here, (j,k),a =t, I is the summation index comprising 
all neighboring sites j and k and the two bosonic species 
't' and J is the single-particle tunneling rate, and 
blj (baj) is the creation (annihilation) operator for a 
boson of type a on lattice site j. The species-independent 
on-site interaction energy is denoted by U, whereas h a j 
are the atomic number operators, 5 a ^> is the Kronecker 
symbol, and \i is the chemical potential. The external 
harmonic confinement with trapping frequency Co>t rap /27r 



^trap 



^mcj t 2 rap (jai at ) 2 , where m is the 



is described by V- 

mass of 8T Rb, and a\ at is the lattice constant. 

In the limit of large interactions (U ^> J) and unity 
filling, it is possible to treat Hbh in second order perturb- 
ation theory [T2HT5] , mapping it to the isotropic spin- 
1/2 Heisenberg Hamiltonian [Eq. (la)], by constructing 
effective spin- 1/2 operators 



(5a) 
(5b) 
(5c) 



and by introducing the exchange coupling J ex = 4J 2 jJJ ' . 

In principle, the parabolic confinement Vi makes su- 
perexchange site-dependent, 



J< 



4J 2 /7 



(6) 



B. Post-selecting samples 

In order to visualize the role of our post-selection 
method, we plotted the probability distributions of 
empty sites obtained from negative images including 
chains with more than one empty site (Fig. SI). We ob- 
serve a clear decrease of the signal contrast as the max- 
imum number of empty sites becomes larger. Although 
this is partially due to the temperature (which we discuss 
in the main text for positive images), the main reason for 
the decrease is that we cannot distinguish between spin 
impurities and excitations within the negative image. 

The post-selection relies on the assumption that there 
is only one spin impurity in the system. This assumption 
is not valid in the following two cases. First, the spin- flip 
on the initial site may fail due to the finite efficiency of 
the multi-site addressing scheme. In this case, if there is 
only one thermal excitation in the region of interest of the 
system, and we count this excitation as the spin impurity. 
Second, if there is a hole at the addressed site, no spin 
impurity is introduced into the system and the hole is 
mistakenly detected as the impurity. A hole is moving 
much faster than the impurity and after averaging over 
samples of different system size its quantum interference 
disappears over the relatively longer time scales on which 
spin- wave dynamics occurs. In both cases, we obtain 
distributions of thermal excitations, which have no clear 
interference-like structure. What we obtained after post- 
selecting samples in our experiment is the sum of the 
position distribution of the spin impurity in the lower- 
temperature systems and the distribution of the thermal 
excitation. This leads to a decrease in contrast, which 
explains the small differences between the experimental 
data and the model in Fig. [4] 

Furthermore, this mechanism for the reduction in con- 
trast is why we claim in the manuscript that the post- 
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However, given the relatively small system size (12 sites) 
and low trapping frequency o; t rap/27r = 85 Hz, the vari- 
ation of J ex across the system does not exceed 2% at 
V x = 10 E r . We can therefore safely consider the super- 
exchange coupling constant, with J ex ,j ~ Jex = 4J 2 //7. 



Figure SI. Post-selecting samples. Shown are position 
distribution of the single-spin impurity for J/U = 0.053 and 
a hold time of 60 ms, after post-selecting samples with one 
(a), less than three (b), four (c), and five (d) empty sites in 
the chain. 
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selected data should actually correspond to an even lower 
temperature (and thus higher contrast) than indicated in 
Fig.[3]i. This remaining ambiguity of the spin-impurity 
position within the post-selected data could in principle 
be resolved by using direct spin-selective single-atom de- 
tection, which would allow for simultaneous imaging of 
both bath and impurity atoms. This, however, is exper- 
imentally challenging. 
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C. Temperature effects 



We investigated the influence of temperature on the 
spin impurity dynamics for J/U = 0.23 (V x = 5E r ), 
in addition to J/U = 0.053 (V x = 10 E r ) which was dis- 
cussed in the main text. For this purpose, we numerically 
simulated the position distributions of the spin impurity 
after an evolution time of 4 ms for different temperatures 
(Fig.|S2fc). The simulation shows that even close to the 
transition point, the propagation velocity as a function 
of temperature remains constant whereas the contrast 
decreases (Fig. S2 j and c), similar to the strongly inter- 
acting limit. 
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Figure S2. Effects of temperature on the dynamics. 

(a) Position distribution of the spin impurity for J/U = 0.23 
and a hold time of 4 ms calculated with t-DMRG for different 
temperatures. Temperature dependence of the contrast (b) 
and the spreading velocity (c) extracted from the dataset of 
(a). 



